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3 Recent research has shown that distinct physical regions in the Venusian 

4 induced magnetosphere are recognizable from the variations of strength and 

5 of wave/fluctuation activity of the magnetic field. In this paper the statis- 

6 tical properties of magnetic fluctuations are investigated in the Venusian mag- 

7 netosheath, terminator, and wake regions. The latter two regions were not 

8 visited by previous missions. We found 1/f fluctuations in the magnetosheath, 

9 large-scale structures near the terminator and more developed turbulence 

10 further downstream in the wake. Location independent short-tailed non-Gaussian 

11 statistics was observed. 
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1. Introduction 

12 In the absence of an intrinsic magnetic field, the structured plasma environment around 

13 Venus is formed due to the direct interaction between the solar wind and the ionosphere. 

14 The ionosphere acts as an obstacle to the supersonic solar wind flow carrying a magnetic 

15 field. The interplanetary magnetic field is draped around the planetary obstacle to form an 

16 induced magnetosphere (Zhang et al., 2007) which consists of the magnetic barrier on the 

17 dayside and the magnetotail. The characteristic scales of the induced magnetosphere and 

18 the magnetotail are different from the spatial dimensions in the Earth's magnetosphere. 

19 The results of Venus Express (VEX) provide a distance of the Venusian subsolar bow shock 

20 from the surface of the planet of about 1900 km, while the terminator bow shock location 

21 is 2.14 Rv {Rv — 6051 km) at solar minimum. The altitude of the Venusian induced 

22 magnetopause is 300 km at the subsolar point and 1013 km at the terminator (Zhang et 

23 al., 2007). For comparison, the width of Earth's magnetosheath is about 25000 km at 

24 the subsolar point. Nevertheless, particularly in case of quasi-parallel shock geometries, 

25 waves and fluctuations occur equally in both Earth's and Venus' magnetosheaths (Luh- 

26 mann et al., 1983, 1986). Due to the rather limited width of the Venusian magnetosheath, 

27 however, the Reynolds number characterizing turbulent flows is presumably small. As a 

28 consequence, the inertial scaling range (if it exists), corresponding to the scales with pre- 

29 vailing nonlinear multi-scale interactions in turbulence, is narrow. Keeping in mind that 

30 both the proton inertial length and the proton Larmor radius are of the order of 100 

31 km, the magnetohydrodynamic (MHD) spatial scales in the Venusian magnetosheath are 

32 limited to one or two decades in wave-number space. For example, the Larmor radius is 
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33 about 160 km for a proton temperature of 1 keV and B=20 nT. The wider scaling region 

34 corresponds to the terminator shock location where the distance between boundaries in- 

35 creases. The terminator and, further downstream, the night-side region is of particular 

36 interest, where plasma instabilities, vortices and wake turbulence can develop due to the 

37 planetary obstacle. The near-polar orbit of VEX with a periapsis altitude of 250-350 km 

38 allows for the first time observations at terminator and mid-magnetotail regions (Zhang 

39 et al., 2006) and thus to compare fluctuation statistics within both magnetosheath and 

40 wake regions. These two important regions were not covered by the previous missions, 

41 e.g. Pioneer Venus Orbiter (Russell, 1992). Using 1-s resolution magnetic field data from 

42 Venus Express during one crossing we provide evidence of changing fluctuation/turbulence 

43 statistical characteristics along the VEX trajectory from the dayside magnetosheath to 

44 the wake (Figure Ic). 

2. Events on May 17, 2006 

45 The magnetic field strength Bt (black line) and the Bx component (grey line) on May 

46 17, 2006 is shown in Figure la. Cartoons of the VEX trajectory (thick black lines) are 

47 shown in VSE (Figure lb) and VSO (Figure Ic) coordinate systems, respectively. In 

48 VSE, the X coordinate is towards the Sun, Z is in the plane of convection electric field 
i9 E — VswxB, {Vsw is the solar wind speed and B is the IMF), and Y completes the 

50 orthogonal coordinate system. The plane of the plasma sheet is ahgned with E (e.g. 

51 Barabash et al., 2007). VEX is crossing the bow shock (thick gray line, quasi-parallel 

52 shock in this case), the magnetosheath, the magnetopause (dashed black line), and the 

53 wake region (approximately the optical shadow - shaded region. Figure Ic). The events on 
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54 this day were initially analyzed by Zhang et al. (2007), who demonstrated that a crossing 

55 of each physical region is recognizable from the variations of strength and variation of the 

56 wave/fluctuation activity of the magnetic field. The magnetic field strength is strongly 

57 fluctuating and its value increases up to ~ 40 nT within the magnetosheath (interval 

58 A in Figure la), after which Bt decreases almost to zero (closest approach, VEX is 

59 below the induced magnetopause) . Magnetic field strength and fluctuations increase again 

60 during interval B (Figures la, b) and this region is defined as the wake. Unlike the 

61 magnetosheath spatial extent, the wake location is more difficult to determine. This 

62 is because, further downstream, VEX is again crossing the magnetopause (at ~ -2 Ry 

63 in Figure Ic), and returns to the magnetosheath and the solar wind (after 03:30 UT 

64 in Figure la). Moreover, the induced magnetopause boundary can be more unstable 

65 downstream. A statistical study of the induced magnetopause boundary crossings (Zhang 

66 et al., 2007) shows that there is a significant uncertainty in the boundary location between 

67 the terminator and X = -2 {Ry)- This indicates that VEX can be located below or over 

68 the flapping magnetopause in an unpredictable manner even during one crossing. In what 

69 follows we refer, for simplicity, to fluctuations downstream of the terminator location, but 

70 before re-entering the solar wind, as wake fluctuations/turbulence. 

2.1. Spectral scalings 

71 We determined the occurrence times of wake fluctuations through spectral analysis. The 

72 spectral scaling indices a, corresponding to magnetic field strength Bt, were estimated 

73 using a wavelet method (Voros et al., 2004) for data intervals of increasing lengths ^2 — ^1, 

74 starting at ti = 01:38 UT (closest approach time). Figure 2 shows that, for t2 extending 
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75 to 0254 UT, two distinct scaling regions appear, exhibiting a sudden transition at about 

76 0208 UT. The robustness of scahng index determinations were further examined for two 

77 sub-intervals before and after the sudden transition at 0208 UT. Optimized data intervals 

78 with smallest errors in a were found. Following this procedure, we identified two distinct 

79 time intervals with well defined a; one between 0139 and 0208 UT and the second between 

80 0209 and 0239 UT. The corresponding time ticks on the orbit are shown in Figure Ic. 

81 The corresponding values of spectral indices during these intervals are ai — 2.45 ± 0.09 

82 and q;2 = 1.65 ± 0.08, respectively (Figure 3). The spectral index in the magnetosheath, 
a — 0.8 ± 0.15, was obtained similarly, for a shorter interval between 0116 and 0131 UT 

84 (filled circles in Figures lb,c). The power spectral density plots in Figure 3 were computed 

85 from the whole magnetosheath interval {A) and two subintervals in the wake {Bi and S2, 

86 see Figure la), both having the length of A. The time intervals A, Bi and B2 are also 

87 indicated by thick black lines along the orbit in Figures lb,c. Bi and B2 are taken from 

88 the two longer, optimized wake intervals, for which the spectral indices ai and a2 were 

89 found previously. The corresponding time ticks {Bi. 0141 - 0156 UT and B2. 0211 - 0226 

90 UT) are shown on the orbit in Figures lb,c. 

91 Figure 3 shows that over the frequency range of 0.04 - 0.5 Hz magnetosheath fluctuations 

92 (thick gray line) exhibit larger power than wake fluctuations. The spectral index a is 

93 close to 1, which indicates the presence of 1// noise. We note that the small peaks 

94 near 0.1 Hz might be associated with mirror mode structures, which occur in planetary 

95 magnetosheaths due to temperature anisotropics (Volwerk et al., 2008, GRL, in press). 
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96 The methods and the data intervals in this paper, however, were not designed for detection 

97 of wavy structures, but for spectral scalings. 

98 The wake fluctuations closer to the terminator (interval Bi , thin black line in Figure 3) 

99 posses significant power only at low frequencies, then the power rapidly decreases following 

100 ~ /~^'^^ (approximately as the thick dashed line shows). Due to the shortness of the time 

101 series the peaks over the frequencies 0.04 - 0.1 Hz are not statistically significant, but 

102 assuming bulk plasma velocity of 200 km/s near the terminator (Spreiter and Stahara, 

103 1992), the largest power is contained in structures of 2000 - 5000 km size. 

104 Numerical simulations indicate that the Kelvin- Helmholtz instability can occur at the 

105 terminator ionopause of Venus (Terada et al., 2002), capable of producing structures over 

106 1000 km and turbulence (Amerstorfer et al.,2007; Biernat et al., 2007). In fact, the 

107 initial VEX observation detected these waves (Balikhin et al., 2008, submitted to GRL). 

108 Our event study shows that MHD turbulence near terminator is not developed because 

109 of the rapid decrease of spectral power toward higher frequences. Further downstream, 
no after 0209 UT (Figure 1) the spectral scahng index is a2 — 1.65 ± 0.08 (thick black line 

111 in Figure 3), very close to 5/3 or 3/2, which are the scaling index values expected for 

112 hydrodynamic or MHD inertial range turbulence. Our spectral analysis shows that the 

113 transition between regions with indices ai and a2 is almost immediate, around 02:09 UT 

114 (see also the waveforms in Figure la). We can speculate, therefore, that the spacecraft 

115 suddenly enters into a wake region with more developed turbulence rather than gradually 

116 developing from the large-scale fluctuations observed near the terminator. The occurrence 

117 of turbulence during the interval B2 might also be related to crossings of the plasma sheet. 
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118 In fact, during this interval Bx is changing sign several times and Bt is decreasing locally 

119 (Figure la), which could be interpreted as a signature of multiple plasma sheet crossings. 

120 Figures lb,c show, however, that the interval B2 is outside of the central part of the wake, 

121 that is, outside of the region where the £^-field ahgned plasma sheet could be situated at 

122 \YVSE\ < O.bRv (e.g. Barabash et al., 2007). 

2.2. Comparison to Gaussian distributions 

123 Higher order statistics is needed to fully describe the nature of nonhnear fluctuations. In 

124 turbulent, boundary influenced non-homogeneous plasma flows, the shape of the probabil- 

125 ity density functions (PDFs) is scale-dependent, non-symmetric and peaked with long tails 

126 (e.g. Voros et al., 2007a). Because of the shortness of our time series the shapes of PDFs 

127 or their scale dependency cannot be evaluated. Instead, we construct Gaussian (normal) 

128 probability plots (GPPs) to decide whether or not our data are Gaussianly distributed. 

129 GPPs are plots of cumulative probabilities versus sorted data values with distorted y axis 

130 labels. The latter ensures a linear dependency for the theoretical Gaussian distributions 

131 in GPPs (e.g. D'Agostino et al., 1990). We reconstructed GPPs using the magnetic field 

132 strength Bt as well as two-point differences defined through 5B = Bxit + r) — Bxit), 

133 for T = 2... 30 s. All these data sets taken from the magnetosheath and wake regions ex- 

134 hibit non-Gaussian statistics. A representative behavior for r = 3 s is depicted in Figure 

135 4. Points are experimental data while the solid lines correspond to theoretical Gaussian 

136 distribution. Near-terminator wake data (interval Bl in Figure la) are the closest to a 

137 Gaussian distribution. Nevertheless, all data sets exhibit nonlinearity, mainly the first few 

138 (negative data values, above the Gaussian line) and the last few (below the Gaussian line) 
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139 experimental points. This is a characteristic behavior for non-Gaussian PDFs with short 

140 tails. Long-tailed data would deviate from the linear relationship in just the opposite 

141 way: the first few points would appear below and the last few points above the Gaussian 

142 solid line (e.g. D'Agostino et al., 1990). The absence of long-tailed PDFs in case of the 

143 turbulent wake (interval B2, Figure la) with scaling index a2 — 1.65 ± 0.08 is surprising, 

144 because the statistics of developed turbulence is usually described in terms of long-tailed 

145 PDFs. The tails of PDFs are built-up from rare but more energetic events, therefore, the 

146 observed short tails can just indicate that our time series are short and not representa- 

147 five enough to observe the expected long-tailed PDF. This assumption can be checked 

148 by investigation of similar crossings and multiple representations of wake turbulence near 

149 Venus. 

3. Conclusions 

150 In this paper the unique data from the VEX spacecraft were used to compare magnetic 

151 fluctuation statistics from Venusian magnetosheath, terminator and mid-magnetotail re- 

152 gions. For simplicity, we denote the latter two regions Venusian wake, never visited during 

153 previous missions. To identify spectral scaling ranges and indices, we used a wavelet tech- 

154 nique, successfully applied for studying the continuous spectra in the Earth's plasma sheet 

155 turbulence (Voros et al., 2004). Our data intervals with the smallest errors in spectral 

156 scaling index corresponded to different physical regions in real space. 

157 We found 1/f non-Gaussian noise in the magnetosheath as signature of the presence 

158 of independent driving sources (see e.g. Voros et al., 2007b). The quasi parallel bow 
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159 shock source for mirror mode waves, the instabihties and gradients near the magnetic field 

160 draping region - can probably contribute to the day side fluctuations in the magnetosheath. 

161 Near the terminator and further downstream in the Venusian wake, large-scale non- 
162 Gaussian structures are present, which can be partially explained in terms of the Kelvin- 

163 Helmholtz instability. As VEX continued further downstream into the wake, more de- 

164 veloped turbulence was detected with spectral scaling index typical for hydrodynamic or 

165 MHD turbulence in the inertial range of scales. The observed non-Gaussian short-tailed 

166 distributions can be caused by the shortness of our data sets. This calls for studying 

167 PDFs from multiple crossings and realizations of turbulent processes. 
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Figure 1. (a) Magnetic field strength Bt (black line) and the Bx component (grey hne) 
on May 17, 2006; the capital letters A, Bl and B2 indicate intervals of equal lengths (black 
lines in each subplot) within the magnetosheath (A) and wake (B); (b) VEX orbit (black 
line) in the VSE coordinate system of the convection electric field E; time ticks on the 
orbit indicate the analyzed intervals; the optical shadow is indicated by the dashed circle, 
the continuous circle shows the terminator bow shock location; (c) VEX trajectory (solid 
black line) in VSO coordinate system; bow shock (solid gray line); induced magnetopause 
(dashed line); optical shadow (shaded region); time ticks on the orbit correspond to the 
analyzed intervals. 
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Figure 3. Power spectral densities for intervals of equal length in the magnetosheath 
(A), and in the wake (Bl and B2); the numerical values of spectral indices were obtained 
during data intervals with the smallest errors. 
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Figure 4. Gaussian (normal) probability plots for data intervals in the magnetosheath 
and wake. The solid lines correspond to the theoretical Gaussian, points correspond to 
the experimental data. 
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